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ABSTRACT 
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An  Earth  impact  with  an  asteroid  has  potentially  devastating  consequences.  In 
order  to  avoid  global  destruction,  astronomers  have  been  attempting  to  map  out  the  orbits 
of  all  possibly  hazardous  Near-Earth  Asteroids  (NEAs)  for  the  past  two  decades. 

Although  there  are  a  number  of  astronomers  devoted  to  discovering  new  NEAs,  there  are 
very  few  who  are  tracking  these  discoveries.  Therefore  many  recently  identified  NEAs 
are  routinely  lost. 

This  Trident  Scholar  Project  had  two  main  objectives.  The  scientific  goal  of  this 
project  focused  upon  tracking  several  NEAs  in  order  to  determine  their  orbits.  In 
addition,  lightcurves  were  observed  for  a  subset  of  these  NEAs  in  order  to  find  their 
rotational  periods  and  their  general  shapes.  This  was  done  by  taking  observations  of  each 
asteroid  for  at  least  one  full  night  and  then  plotting  light  intensity  versus  time. 

The  second  objective  was  to  refurbish  and  upgrade  the  observatory  at  the  U.S. 
Naval  Academy.  Astronomical  research  at  the  U.S.  Naval  Academy  hinges  on  the 
successful  operation  of  a  20"  Cassegrain  telescope  with  a  Photometries  1024x1024 
Charge  Coupled  Device  (camera),  as  well  as  several  computers  to  control  these  devices 
and  process  the  data.  As  a  result  of  this  Trident  Scholar  Project,  the  U.S.  Naval  Academy 
now  boasts  a  fully  functioning,  state-of-the-art  observatory,  including  newly  acquired  and 
developed  software  for  tracking  NEAs  and  analyzing  astronomical  observations  in 
general. 

This  Trident  Scholar  Project  added  to  the  body  of  knowledge  of  NEAs  by 
determining  the  orbits  of  eleven  asteroids,  as  well  as  the  rotational  periods  and  rough 
shapes  of  four  asteroids.  In  addition,  this  Trident  Scholar  Project  will  benefit  both 
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midshipmen  and  faculty  who  wish  to  either  continue  tracking  NEAs  or  begin  a  new  field 
of  astronomical  research  using  the  now  completely  established  observatory  at  the  U.S. 
Naval  Academy. 

KEYWORDS:  NEA,  Near-Earth  Asteroid,  Asteroid,  Telescope,  Charge-Coupled 


Device,  CCD,  Orbit,  Lightcurve,  Astronomy 


ACKNOWLEDGEMENTS 


3  ■ 


I  would  like  to  thank  Assoc.  Prof.  Katz-Stone  for  all  of  her  help,  guidance,  and 
long  nights  throughout  this  project.  I  would  also  like  to  thank  Dr.  Jeffery  Larsen  at 
Spacewatch  for  his  help  in  understanding  how  to  do  this  project.  I  would  like  to  thank 
Prof.  Albert  for  the  advice  and  help  she  has  given  me  to  get  the  observatory  up  and 
running.  I  would  like  to  thank  Jure  Skvarc  for  his  assistance  with  Fitsblink.  I  would  also 
like  to  thank  Dave  Moore,  Dave  Rector,  and  Jeff  Walbert  for  helping  me  take  care  of  the 
telescope.  I  would  like  to  thank  Ian  Huss  and  Mark  Kelly  for  their  help  with  fixing  the 
telescope’s  many  problems.  I  would  finally  like  to  thank  the  Trident  Scholar  Committee, 
Dean  Miller,  and  the  United  States  Naval  Academy  Physics  Department  for  allowing  me 
to  do  this  project. 


TABLE  OF  CONTENTS 


4 


ABSTRACT _  1 

ACKNOWLEDGEMENTS _ 3 

1  INTRODUCTION _ 5 

2  OBSERVATIONS _ 8 

2.1  Equipment  8 

2.2  Determining  General  Instrument  Parameters  12 

2.3  A  Typical  Night’s  Observation  15 

2.4  Calibration  and  Data  Reduction  1 7 

2.5  Identifying  NEAs  18 

3  ORBITS _ 20 

3.1  Gauss’  Method  22 

3.2  Determination  of  Orbital  Elements  27 

3.3  Perturbed  Motion  29 

4  LIGHTCURVES _ 30 

5  RESULTS  AND  DISCUSSION _ 32 

5.1  Cooperation  Process  32 

5.2  Using  an  Astronomical  Calibrator  35 

5.3  Orbital  Results  and  Discussion  36 

5.4  Lightcurve  Results  and  Discussion  49 

6  CONCLUSION _ 55 

7  FUTURE  WORK _  56 

REFERENCES _ 57 

APPENDIX  A:  LAPLACE’S  METHOD  58 


1.  Introduction 


5 


In  1993,  Drs.  Carolyn  and  Eugene  Shoemaker  and  amateur  astronomer  David 
Levy  made  an  astonishing  discovery.  They  found  a  fragmented  comet,  later  named 
Shoemaker-Levy  9,  caught  in  Jupiter’s  gravitational  field.  Over  the  course  of  about  a 
year,  astronomers  around  the  world  watched  as  the  comet  further  split  apart  into  more 
than  twenty  fragments  and  then  as  each  fragment  struck  Jupiter.  The  impacts  ranged 
from  barely  discernable  to  the  equivalent  of  about  six  million  megatons  of  TNT.  If  comet 
Shoemaker-Levy  9  had  hit  Earth,  it  is  unlikely  there  would  be  any  human  life  remaining 
(e.g.,  Zelik  &  Gregory  1998;  Boslough  &  Crawford  1995  and  references  therein). 

In  fact,  it  is  likely  that  such  an  impact  destroyed  the  dinosaurs.  Luis  and  Walter 
Alvarez  studied  an  impact  crater  on  the  Yucatan  Peninsula  (Alvarez  et  al.  1980).  They 
found  a  layer  of  iridium  in  the  Cretaceous-Tertiary  boundary,  which  corresponds  to  the 
epoch  of  dinosaur  extinction.  Iridium  is  not  normally  found  on  the  Earth’s  crust  because 
it  is  believed  to  have  sunk  to  the  center  of  the  Earth  during  its  formation;  on  the  other 
hand,  asteroids  are  often  rich  in  iridium.  Therefore,  the  Alvarezes  believe  that  the 
iridium  must  have  come  from  a  large  asteroid  or  comet  (Alvarez  et  al.  1980;  Alvarez 
1990),  and  this  impact  led  to  the  extinction  of  the  dinosaurs.  Studies  of  these  two 
powerful  impacts  have  awoken  an  increasing  sense  of  urgency  in  the  scientific 
community  and  in  the  U.S.  Congress.  Accordingly,  a  new  goal  has  been  set  to  map  the 
orbit  of  all  2000  potentially  devastating  near-earth  objects. 

The  size  of  the  asteroid  determines  how  destructive  it  will  be  upon  impact.  An 
asteroid  that  is  10  km  in  diameter  will  cause  earthquakes  of  magnitude  12  on  the  Richter 
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Scale1  with  amplitudes  of  200  m.  The  energy  from  this  sort  of  impact  will  raise  the 
atmosphere’s  temperature  to  about  500°  F.  Large  chunks  of  rock  would  be  thrown  up  in 
the  air.  Many  would  escape  the  pull  of  the  Earth’s  gravity.  The  rest  would  fall  back  to 
Earth,  crashing  into  the  entire  planet.  An  asteroid  of  about  1  km  in  diameter  could 
destroy  a  continental  area  directly,  and  it  would  cause  a  nuclear  winter  due  to  the  dust 
thrown  up  into  the  atmosphere.  The  dust  would  block  the  Sun’s  radiation,  cooling  the 
Earth  to  temperatures  far  too  cold  for  humans.  An  asteroid  of  about  0.1  km  would  not 
annihilate  the  entire  world  population,  but  it  would  cause  a  great  amount  of  damage.  If  it 
hits  land,  an  asteroid  of  this  size  could  easily  destroy  an  area  the  size  of  a  state.  If  a  0.1 
km  asteroid  hits  water,  which  is  very  likely  because  the  surface  area  of  the  Earth  is  3/4 
water,  it  would  create  a  tsunami  that  would  drown  an  entire  coastline.  Finally,  a  10  m 
asteroid  poses  no  threat  to  humanity  as  a  whole.  However,  the  crater,  debris,  and  fires 
caused  by  its  impact  would  still  cause  much  damage  (see,  e.g.  Amy  1998;  Toon  et  al. 
1994). 

In  order  to  focus  upon  preventing  global  destmction,  it  is  imperative  to  find  and 
track  all  Near-Earth  Asteroids  (NEAs)  that  are  at  least  1  km  in  diameter.  This  world¬ 
wide  effort  is  coordinated  by  the  Minor  Planet  Center2  which  collects  data  on  each  newly 
discovered  asteroid  and  disseminates  the  necessary  information  via  the  world-wide  web 


1  The  Richter  Scale  is  a  logarithmic  scale,  devised  in  1935  to  identify  the  magnitude  of  an  earthquake.  The 
value  calculated  reflects  the  amount  of  energy  released.  The  Richter  scale  ranges  from  0  to  8.9,  the  highest 
recorded  quake.  Because  it  is  a  logarithmic  scale,  an  earthquake  of  6.0  is  10  times  more  severe  than  one  of 
5.0. 

2  The  Minor  Planet  Center  (MPC)  operates  at  the  Smithsonian  Astrophysical  Observatory,  under  the 
auspices  of  Commission  20  of  the  International  Astronomical  Union  (IAU)  and  is  a  nonprofit  organization, 
with  principal  funding  coming  from  subscriptions  to  the  various  services  offered  by  the  Center.  The  MPC  is 
responsible  for  the  efficient  collection,  computation,  checking  and  dissemination  of  astrometric 
observations  and  orbits  for  minor  planets  and  comets,  via  the  Minor  Planet  Circulars  (issued  monthly)  and 
the  Minor  Planet  Electronic  Circulars. 
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to  other  astronomers  who  do  the  follow-up  tracking  of  these  newly  discovered  asteroids. 
Only  a  handful  of  observatories3  are  dedicated  to  discovering  new  asteroids  and  even 
fewer  professional  observatories  track  asteroids.  In  fact,  the  majority  of  NEAs 
discovered  are  lost.  The  major  goals  of  this  Trident  Project  were  to  develop  a  program 
here  at  USNA  to  track  NEAs  and  to  observe  several  NEAs  in  order  to  determine  their 
orbits,  rotational  periods,  and  basic  shapes. 


3  In  the  United  States  this  effort  is  led  by  Spacewatch  in  Arizona,  LINEAR  in  New  Mexico  and  NEAT  in 
Hawaii. 


2.  Observations 
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2. 1  Equipment 

The  telescope  used  for  this  project  is  a  Cassegrain  reflecting  telescope 
manufactured  by  DFM  Engineering,  Inc.4  and  mounted  on  the  roof  of  Michelson  Hall 
(Figure  1).  The  telescope  has  a  20"  parabolic  primary  mirror,  which  collects  the  light 
from  astronomical  sources  and  focuses  it  on  a  convex  secondary  mirror  in  the  middle  of 
the  telescope.  This  secondary  mirror  then  reflects  the  light  through  a  hole  in  the  center  of 
the  parabolic  mirror  and  focuses  it  on  the  detector  (e.g.,  Shu  1982).  In  addition,  an  IBM 
computer  controls  the  movements  of  the  telescope  and  dome.  In  particular,  the  dome  and 
telescope  must  spin  to  compensate  precisely  for  the  Earth’s  rotation  rate.  This  allows 
observations  of  a  single  object  for  up  to  an  hour  without  resetting  the  telescope's  pointing 
position. 

A  good  astronomical  detector  is  a  device  that  produces  a  predictable  response  to 
the  light  gathered  and  focused  by  the  telescope.  Early  astronomers  used  the  human  eye 
as  their  detector.  This  was  later  replaced  by  photographic  film  and  photomultipliers. 
Today,  modem  astronomers  use  charge  coupled  devices  (CCDs),  which  were  first 
developed  in  the  early  1980's.  While  their  costs  may  be  prohibitive  for  some  observers, 
their  advantages  over  other  detectors  are  crucial  for  many  projects.  First,  CCDs,  unlike 
photographic  film,  have  a  linear  response  rate5  over  a  large  range  of  intensities.  Second, 
CCDs  are  very  sensitive,  detecting  about  50%  of  the  incident  photons.  In  addition,  they 

4  DFM  Engineering,  Inc.  Longmont,  Colorado 

5  The  CCD’s  response  is  directly  proportional  to  the  number  of  incident  photons. 


Fig.  1.—  This  DFM  telescope  was  used  for  all  the  observations  taken  at  USNA.  The  small  telescope  is  the 
finder  scope. 


Fig.  2.  —  Left:  A  close-up  of  the  Photometries  CCD  used  to  image  all  the  observations  taken  at  USNA. 
Right:  The  interior  of  the  CCD.  The  small  square  near  the  center  of  the  image  is  the  location  of  the  actual 
CCD  chip. 
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are  high  resolution  detectors  with  micron-sized  pixels.  Another  advantage  is  that,  with 
a  CCD,  the  contamination  caused  by  nearby  stars  or  cosmic  rays  can  be  removed  from 
the  image.  A  final  advantage  is  that  the  readout  of  a  CCD  is  already  digitized  and  easily 
fed  directly  into  a  computer,  which  is  used  to  calibrate  and  analyze  the  data. 

The  detector  used  for  this  project  is  a  CCD  manufactured  by  Photometries6 
(Figure  2).  Like  all  CCDs,  it  is  made  up  of  photosensors,  which  can  be  imagined  as  little 
squares  on  a  chessboard  but  so  small  that  millions  can  fit  in  a  few  square  centimeters. 

Our  CCD  has  1024  by  1024  photosensors.  Each  photosensor  is  made  of  a  light  sensitive 
material  which,  when  struck  by  a  photon,  releases  an  electron.  Each  pixel  is  held  at  a 
positive  potential  so  that  the  electrons  are  trapped  in  the  potential  well  during  the 
exposure.  The  CCD  is  controlled  by  a  Micron  computer,  which  is  also  used  to  archive 
and  calibrate  the  data. 

After  the  exposure  is  completed,  the  computer  reads  the  potential  off  each  pixel 
one  at  a  time.  The  CCD  does  this  with  a  horizontal  register,  which  has  as  many  pixels  as 
a  line  in  the  array  (1024)  but  is  covered  so  as  not  to  be  exposed  to  the  light.  In  order  to 
read  the  potential  off  of  the  CCD,  vertical  clocks  transfer  each  line  in  a  CCD  to  the  one 
directly  below  it,  and  horizontal  clocks  transfer  the  pixel’s  electric  charge  laterally  down 
the  horizontal  register  to  the  register’s  extremity  for  the  output,  as  shown  in  Figure  3. 

This  procedure  is  repeated  until  all  of  the  pixels  have  been  read  (e.g.,  Martinez  &  Klotz 
1999).  Once  this  is  done,  the  computer  converts  the  potential  of  each  pixel  into  counts 
and  assembles  a  picture. 


6  Photometries,  Roper  Scientific  Inc.  Tucson,  Arizona. 


Fig.  3.  —  This  illustrates  how  a  single  pixel  is  read  off  the  CCD.  This  process  repeats  until  all  pixels 
have  been  read. 
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The  CCD  is  exposed  to  light 
and  a  charge  pattern  forms. 


The  charge  is  shifted  one  column. 
The  first  column  is  shifted 
into  the  serial  register. 


The  charge  in  the  serial 
register  is  shifted  up.  The 
first  pixel  is  shifted  into 
the  output  node  where  it 
is  subsequently  collected  for 
processing. 
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Unfortunately,  the  CCD  is  not  only  sensitive  to  the  light  from  the  secondary 
mirror,  but  also  to  thermal  radiation  produced  by  the  electronics  of  the  CCD  itself,  known 
as  the  “dark”  current.  To  reduce  the  dark  current,  antifreeze  is  used  to  chill  the  CCD  to 
approximately  -40°C.  This  is  further  complicated  because  the  dark  current  increases  the 
longer  a  pixel  is  exposed  to  the  thermal  radiation.  So  as  each  row  of  the  CCD  is 
transferred  to  the  horizontal  register,  it  is  still  exposed  to  the  dark  current.  The  first  row 
read  has  a  smaller  dark  current  than  the  last  row  read.  Therefore,  all  images  have  a 
gradient  across  them. 

Finally,  the  sensitivity  of  each  pixel  varies.  The  response  factor  for  the  pixels  can 
be  affected  by  the  gain  irregularities  of  the  CCD’s  pixels  and  from  problems  with  optical 
transmission,  such  as  vignetting7  and  the  presence  of  dust  in  the  light  path  (Martinez  and 
Klotz  1999). 

2,2  Determining  General  Instrument  Parameters^ 

Before  observing  any  asteroids,  some  basic  information  about  the  equipment  was 
needed.  Between  the  telescope  and  detector  is  an  eight-filter  wheel.  Each  filter  allows 
only  one  waveband  of  light  through. 

The  focal  point  of  the  optical  system  is  different  for  each  filter.  The  focus  for 
each  filter  is  determined  by  trial  and  error.  For  each  filter,  a  series  of  pictures  of  a  star 
field  were  observed  at  different  focuses.  The  image  in  which  the  stars  appear  point-like 
rather  than  donut-like  corresponds  to  the  right  focus  for  that  filter.  In  order  to  choose  the 
best  filter  for  observing  these  faint  NEAs,  NGC  7790,  a  known  star  cluster  was  observed 
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using  all  eight  of  the  filters  as  well  as  without  a  filter.  NGC  7790,  as  shown  in  Figure 
4,  is  a  cluster  of  stars  whose  magnitudes,  image,  and  scale  has  been  previously  observed 
and  published  (Odewahn  et  al.  1992).  The  best  images  of  NGC  7790  were  observed  with 
a  Kodak  filter  that  removes  any  light  bluer  than  5900A8.  This  filter  removes  artificial 
light  produced  by  sodium  and  mercury  lamps. 

Next,  the  resolution  and  field  of  view  of  the  instrument  needed  to  be  determined. 
This  was  done  by  observing  the  same  standard  star  field,  NGC  7790.  The  magnitudes  of 
the  stars  and  distances  between  the  stars  on  the  plane  of  the  sky  are  known.  The 
resolution  is  determined  by  taking  a  number  of  measurements  of  the  distance  between 
stars  on  both  the  exposure,  Figure  5,  in  pixels  and  the  published  picture  of  NGC  7790  in 
arcseconds9.  Comparing  these  two  measurements  gives  the  pixel  size  in  arcseconds.  By 
taking  several  measurements,  they  can  be  averaged  to  obtain  a  reliable  value  of  0.885 
arcsec/pixel.  The  diffraction10  limited  resolution  of  a  20"  telescope  is  0.23  arcsec.  Since 
this  is  so  much  smaller  than  the  pixel  size,  the  pixel  size  effectively  becomes  the 
resolution.  The  field  of  view  is  then  just  1024  pixels  x  0.885  arcsec/pixel  =  906  arcsecs. 
This  is  about  half  of  the  angular  size  of  the  moon. 

Finally,  this  same  standard  field  was  used  to  determine  the  exposure  time  required 
to  observe  objects  as  faint  as  a  typical  NEA,  about  18  magnitude  .  To  see  the  faintest 
object  in  NGC  7790,  about  18.5  in  magnitude,  a  one-minute  exposure  was  needed. 

7  Light  loss  or  “fall-off 5  from  the  center  of  the  lens  to  the  outside  edge  due  to  the  lens  being  slightly  smaller 
than  the  array. 

8  lA=10'10m 

9  1  arcsec=  1/3600  degrees 

10  Diffraction  describes  the  bending  of  light  as  it  passes  through  an  aperture. 

1 1  A  unitless  number,  based  on  a  logarithmic  scale,  used  to  describe  the  brightness  of  an  astronomical 
object.  Faint  objects  have  larger  magnitudes  than  bright  objects 
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Fig.  4.  —  This  star  field  (NGC  7790)  was  observed  by  Odewahn  et  al.  (1992).  It  was  used  to  determine 
basic  parameters  for  optical  system. 


mi-m2=2.51ogf2/fi 

where  fj  is  the  flux  (power/area)  of  object  i  and  m*  is  its  magnitude. 


Fig.  5.  —  This  image  of  NGC  7790  was  observed  with  the  USNA  observatory.  Stars  in  Figure  4  and  the 
boxed  region  were  compared  to  determine  basic  parameters  as  explained  in  the  text. 

2.3  A  Typical  Night ’s  Observation 


A  typical  night’s  observation  begins  in  the  late  afternoon  when  the  CCD  and 
chiller  are  turned  on  so  that  the  CCD  has  enough  time  to  be  cooled  to  an  operating 
temperature  of  -40°C.  Then,  just  before  dusk,  flat-field  images  are  taken  which  will  be 
needed  later  to  calibrate  the  data.  Flat-field  images  are  used  to  correct  for  differences  in 
the  response  of  the  photosensors.  To  reduce  the  dark  current  in  these  images,  the 
exposure  time  is  short,  less  than  five  seconds.  To  increase  the  signal-to-noise  ratio  (SNR) 
in  the  final  flat-field  image  and  to  reduce  the  effects  of  the  transient  events  such  as 
cosmic  rays,  twenty  to  forty  flat-field  images  are  taken,  and  the  median  of  these  images  is 
found.  Flat-field  images  should  be  uniformly  illuminated  and  free  of  stars.  Therefore,  as 
the  sun  sets,  the  telescope  is  moved  to  the  west  to  avoid  imaging  stars. 
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About  a  half-hour  before  observing  an  asteroid,  the  telescope’s  alignment  needs 
to  be  calibrated.  The  telescope  is  manually  pointed  at  a  known  object,  and  a  picture  is 
taken.  Usually,  the  image  is  not  centered.  The  telescope  is  manually  adjusted  until  the 
object  is  at  the  center  of  the  CCD.  Now,  the  telescope’s  computer  is  reinitialized.  This  is 
a  difficult  process  because  the  CCD’s  field  of  view  is  too  small  to  identify  constellations 
and  individual  stars  look  alike.  In  order  to  overcome  this  difficulty,  the  telescope  is  first 
pointed  at  a  planet  such  as  Jupiter  or  Saturn,  a  very  distinguishable  object.  Once  the 
planet  is  centered  and  the  computer  is  reinitialized,  the  telescope  is  close  to  being  aligned. 
To  refine  the  alignment,  the  computer  is  used  to  point  the  telescope  to  a  bright  star  such 
as  Vega.  A  short  exposure  is  taken  and  only  Vega  appears  near  the  center  of  the  CCD, 
because  the  alignment  is  close  and  Vega  is  so  bright.  Now,  the  telescope  is  manually 
moved  until  Vega  is  in  the  center  of  the  CCD  and  the  computer  is  reinitialized  one  last 
time.  That  aligns  the  telescope  accurately. 

Finally,  an  NEA  can  be  observed.  The  Minor  Planet  Center  maintains  a  list  of 
newly  discovered  asteroids  on-line.  It  also  provides  software,  which  extrapolates  an 
asteroid’s  ephemeris.12  The  telescope  is  pointed  at  the  asteroid  and  an  exposure  is  taken. 
In  a  single  image,  an  asteroid  looks  exactly  like  a  faint  star.  In  order  to  distinguish  an 
asteroid  from  a  star  at  least  two  images  must  be  taken;  in  principle,  the  object  that  moved 
between  exposures  is  the  asteroid.  The  asteroid  must  move  fifty  pixels  before  another 
exposure  can  be  taken  so  that  the  asteroid  has  moved  far  enough  for  the  human  eye  to 
easily  distinguish  the  change  in  position.  The  time  between  exposures  is  determined 
from  the  estimated  speed  of  the  asteroid.  While  only  two  exposures  are  required  to  detect 


12  An  ephemeris  is  the  estimated  position  of  an  object  as  a  function  of  time  based  on  previous  observations. 
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an  asteroid,  at  least  three  exposures  are  taken  for  each  asteroid  in  case  there  is  a 
problem  with  one  of  the  images  and  to  help  properly  identify  the  asteroid. 

Before  the  end  of  the  night,  two  more  series  of  images,  which  are  used  to  calibrate 
the  data,  must  be  taken.  As  described  in  Section  2.1,  each  photosensor  has  an  initial 
positive  potential.  To  measure  this  pre-charge  potential,  a  series  of  “bias”  images  is 
taken.  These  are  short  (1  ms)  exposures  taken  with  the  shutter  closed.  The  other  series 
of  images  are  used  to  measure  the  dark  current.  The  dark  current  depends  on  both  the 
temperature  of  the  CCD  and  the  exposure  time.  Therefore,  the  series  of  dark  images 
must  have  the  same  exposure  time  that  was  used  during  the  asteroid  observations  and 
with  the  CCD  at  the  same  temperature. 

2.4  Calibration  and  Data  Reduction 

As  explained  in  Sections  2.1  and  2.3,  a  raw  image  is  contaminated  by  the  dark 
current  and  pre-charge  potential  and  has  a  gain,  which  varies  from  pixel  to  pixel. 
Calibration  is  the  process  that  corrects  for  these  known  effects.  In  a  raw  image,  the 
intensity  of  any  pixel  (x,y)  is  given  by 

I(x,y)=b(x,y)+d(x,y,t,T)+i(x,y)r(x,y)  ( 1 ) 

where  I(x,y)  is  the  raw  image,  b(x,y)  is  the  bias  value,  d(x,y,t,T)  is  the  dark  value,  r(x,y) 
is  the  response  factor,  and  i(x,y)  is  the  final  calibrated  image  for  an  exposure  time  t  at 
temperature  T  (e.g.,  Martinez  and  Klotz  1999).  Solving  for  the  calibrated  image,  gives 
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The  bias  value  is  estimated  by  taking  a  series  of  very  short  exposure,  with  the 
shutter  closed.  Since  the  exposure  time  is  short,  the  dark  current  built  up  during  the 
exposure  can  be  ignored.  To  increase  the  SNR  in  the  bias  image,  the  median  of  a  series 
of  ten  bias  images  is  found, 

lb(x,y)=b(x,y).  (3) 

where  lb  is  the  median  bias  image.  A  median  is  taken,  not  a  mean,  because  the  normal 
bias  value  for  each  pixel  is  being  determined.  However,  transient  events  such  as  cosmic 
rays  affect  the  bias  value  of  a  random  pixel.  By  taking  a  median,  this  abnormality  is 
ignored.  If  a  mean  is  taken,  however,  a  cosmic  ray  would  dominate  the  average  bias 
value  of  that  pixel  and  thus  gives  the  incorrect  bias  value. 

The  dark  value  is  found  by  taking  an  image  with  the  shutter  closed  at  the  same 
temperature  and  for  the  same  exposure  time  as  the  image.  As  with  the  bias  images,  the 
median  is  found  from  a  series  of  ten  images.  However  this  median  image  (la)  also 
contains  the  bias  value,  so 

ld(x,y)=b(x,y)+d(x,y,t,T).  (4) 

Finally,  the  response  factor  is  found  by  taking  a  series  of  flat-field  images  (as 
described  in  Section  2.3)  and  finding  their  median.  The  flat-field  image  still  contains  the 
bias  value,  but  has  an  exposure  time  of  one  to  five  seconds  so  that  the  dark  current  can  be 
ignored. 

lf(x,y)=b(x,y)+r(x,y)*C  (5) 

The  normalization  constant,  C  is  the  average  of  the  median  flat-field  image.  Substituting 
equations  (3),  (4),  and  (5)  into  (2)  gives  the  formula  for  the  calibrated  image, 


i(x,y)=^±.k(*-y>.c 

lf(x,y)-lb(x,y) 


(6) 
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2.5  Identifying  NEAs 

Once  the  calibrated  image  is  determined,  the  next  step  is  identifying  the  asteroid 
using  an  analysis  program  known  as  Fitsblink13  running  on  a  Dell  Computer.  First,  all  of 
the  real  sources  in  the  image  are  determined.  The  parameters  must  be  chosen  carefully  so 
that  all  the  real  sources  are  identified  without  picking  up  false  detections.  For  example, 
real  sources  tend  to  be  at  least  two  pixels  in  extent  and  have  at  least  900  counts.  NEAs, 
unlike  stars,  move  rapidly  across  the  sky.  Therefore  the  easiest  way  to  distinguish  an 
asteroid  from  a  faint  star  is  to  observe  its  motion.  This  is  done  by  taking  two  or  more 
images.  As  explained  above,  the  time  separation  of  the  two  images  should  be  such  that 
the  asteroid  moves  fifty  pixels  in  that  time.  Then,  the  two  images  are  registered  on  top  of 
one  another  by  aligning  the  same  bright  stars  found  in  both  images.  Next,  the  two  images 
are  “blinked”,  which  means  the  computer  display  rapidly  switches  between  the  two 
images.  Since  the  images  are  aligned,  the  object  that  moves  must  be  the  asteroid.  Once 
the  asteroid  has  been  detected,  the  stars  in  the  field  need  to  be  matched  to  stars  in  the 
Guide  Star  Catalog14  or  the  U.S.  Naval  Observatory  Catalog15  in  order  to  get  the  position 
and  velocity  (astrometry)  as  well  as  the  magnitude  (photometry)  of  the  asteroid.  Once 


13  Fitsblink  is  a  program  written  by  Jure  Skvarc  that  was  used  to  do  astrometry  and  photometry.  It  aligns 
two  images  over  one  another  and  blinks  them  to  identify  the  asteroid.  It  also  compares  and  matches  an 
observer's  image  to  an  image  from  the  Guide  Star  Catalog  so  that  the  coordinates  of  every  pixel  are 
determined. 

14  The  guide  star  index  is  provided  by  the  Hubble  Space  Telescope  Science  Institute.  It  was  constructed 
from  1477  Schmidt  plates;  it  contains  entries  for  over  18  million  objects  between  8.5  and  15th  magnitude. 

15  The  USNO  SA2.0  guide  star  index  was  put  together  by  the  U.S.  Naval  Observatory.  It  contains  over  50 
million  stars  observed  in  both  the  blue  and  red  Palomar  Optical  Sky  Survey. 
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the  astrometry  is  completed,  a  new  ephemeris  is  calculated  so  that  the  next 
observational  position  is  determined,  and  these  data  will  be  used  in  the  final  calculation 


of  the  asteroid’s  orbit. 


3.  Orbits 
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Near-Earth  Asteroids  (the  focus  of  this  project)  are  split  into  three  categories: 
Amors,  Apollos,  and  Atens.  Amor  asteroids  have  orbits  that  are  inside  the  orbit  of  Mars. 
While  Amors  never  cross  Earth’s  orbit,  they  may  still  travel  close  to  Earth's  orbit.  Apollo 
asteroids  are  sometimes  called  Earth-crossing  asteroids  because  their  aphelion16  is 
beyond  Earth’s  orbit  while  their  perihelion17  is  inside  the  Earth’s  orbit.  Aten  asteroids 
have  an  orbit  with  a  semi-major  axis  close  to  1  AU18,  so  they  spend  most  of  their  time  in 
the  Earth’s  orbital  path  (e.g.,  Fix  1999). 

A  single  observation  of  an  asteroid  consists  of  determining  its  position  in  the 
plane  of  the  sky,  not  its  distance  along  our  line  of  sight,  at  the  time  of  observation.  In 
order  to  determine  its  orbit,  a  minimum  of  three  observations  is  required.  With  three 
observations,  the  asteroid’s  heliocentric  position  (r, ),  velocity  (r,. )  and  acceleration  (r, ) 
as  a  function  of  time  are  estimated.  An  orbit  may  be  calculated  from  these  estimates. 

A  fully  determined  orbit  is  described  mathematically  by  six  orbital  elements:  the 
longitude  of  the  ascending  node19  (Q),  the  angle  of  inclination  (/)  with  the  ecliptic20 
plane,  the  argument  of  perihelion  (co),  the  semi-major  axis  (a),  the  eccentricity  (e),  and 
the  mean  anomaly  (Mo).  The  longitude  of  the  ascending  node  (Q)  is  the  angle  measured 


16  An  object's  farthest  point  from  the  Sun. 

17  An  object's  closest  point  to  the  Sun. 

18  An  Astronomical  Unit,  AU,  is  the  distance  from  Earth  to  the  Sun.  1  AU  ~  93  million  miles 

19  The  ascending  node  is  the  point  in  the  orbit  that  crosses  the  ecliptic  from  South  to  North. 

20  The  ecliptic  is  the  plane  of  Earth’s  orbit  and  is  traced  out  by  the  apparent  path  of  the  Sun  during  the  year. 
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9 1 

eastward  from  the  vernal  equinox  to  the  ascending  node.  The  inclination  angle  (i)  is 
the  angle  between  the  ecliptic  and  the  plane  of  orbit  such  that  0  <  i  <180°.  The  argument 
of  perihelion  (co)  is  the  angle  measured  in  the  direction  of  the  orbital  motion,  from  the 
ascending  node  to  the  direction  of  the  perihelion.  These  relationships  can  be  seen 
graphically  in  Figure  6  (e.g.,  Danby  et  al.  1995).  The  semi-major  axis  (a)  is  half  the 
length  of  the  major  axis  of  an  ellipse.  The  eccentricity  (e)  of  a  closed  orbit  is  a  parameter 
between  0  and  1  such  that  a  circle  has  an  eccentricity  of  0.  The  mean  anomaly  (Mo) 
specifies  the  position  of  the  object  at  a  fiducial  time. 


HORTET 


Fig.  6.— Elements  of  an  orbit. 

There  are  two  classical  techniques  for  determining  these  orbital  elements: 
Laplace's  Method  and  Gauss'  Method.  Laplace's  Method  directly  solves  the  equation  of 
motion  described  by  the  gravitational  interactions  in  the  solar  system.  In  practice,  this 
method  is  difficult  because  it  requires  numerical  differentiation  of  a  minimal  amount  of 
data.  Like  the  Spacewatch  observatory,  this  Trident  Project  uses  Gauss’  Method.  Gauss' 


21  The  vernal  equinox  is  the  ascending  node  on  the  celestial  sphere  where  the  celestial  equator  crosses  the 
ecliptic. 


23 

Method  does  not  use  the  equations  of  motion  at  one  instant  of  time.  Instead,  it  uses  the 
solution  of  the  unperturbed  two-body  problem  to  bridge  the  time  interval  between 
observations.  This  method  requires  Taylor  Series  expansions  as  explained  below  (Taff 
1984).  Following  Danby  (1988),  Gauss’  Method  is  reviewed  below  while  Laplace’s 
Method  reviewed  in  Appendix  A. 


5.7  Gauss'  Method 

Gauss'  Method  was  developed  two  hundred  years  ago  and  since  then  has  been 
continually  improved  upon.  It  involves  taking  three  observations,  bridging  them  with  a 
Taylor  Series  expansion  and  solving  for  the  position  rj  of  the  object  at  the  three  times. 
Once  rj  is  found,  the  orbital  elements  can  be  determined. 

As  shown  in  Figure  7,  the  positions  at  the  three  times,  iq,  Y2,  and  r3  all  lie  in  a 
plane  and  Y2  can  be  expressed  in  terms  of  ri  and  r3, 

r2=Cir1+c3r3  (7) 

where  cj  and  C2  are  constants. 


Fig.  7.-  Geometrical  representation  of  equation  (7). 
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The  position  r(t)  can  be  written  as  a  Taylor  series,  which  can  be  reduced  to  an 
expression  in  terms  of  a  fiducial  position  and  velocity  (ro  and  r0 ), 

r(t)  =  f(t)ro+g(t)fo  (8) 

where  ro  and  r0  are  the  position  and  velocity  vectors  at  time  to,  and  f  and  g  are  series, 


+^<#-03  +■■• 

where 


and  P  =  ^ 


(10) 


Therefore,  the  positions  at  ti  and  t3  are  written  in  terms  of  r2  and  r2  as 


r,  =fir2+glr2 
r3  =f3r2  +g3r2, 


(11) 


and 

r2  =Ci(fir2+g,r2)+c3(f3r2+g3r2).  (12) 

The  next  step  is  finding  ci  and  C3  by  eliminating  r*  and  fj .  Taking  the  cross  products  of 
r2  and  r2  with  equation  (12)  yields 

l  =  c,f1+c3f3  and  0  =  c,g,+c3g3.  (13) 


Solving  equation  (13)  for  ci  and  C2  gives 


c,  = 


§3 


flg3-glf3 


c3  =- 


gl 


flg3"glf3 


(14) 
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The  geometrical  association  between  the  Sun,  Earth,  and  the  asteroid  can  be 
seen  in  Figure  8.  Since  r  is  neither  known  nor  even  observed,  the  observer’s  relation  to 
the  asteroid  must  be  introduced,  p.  From  the  geometry  shown  in  Figure  8, 

R  +  r  =p  =p  p.  (15) 

Earth 


r 

Fig.  8.  —  Geometrical  representation  of  equation  (15). 


Substituting  this  into  equation  (7), 

C,P1P1-P2P2+C3P3P3  =C1R1-R2+C3R3  (16) 

and  substituting  in  equation  (1 1) 

fir2  +  gl*2  =  PlPl  -R1 

r2=P2P2-R2  (17) 

f3r2+g3r2  =P3P3“R3- 

In  order  to  solve  for  X2,  it  is  convenient  to  change  coordinate  systems  from  a 
heliocentric  system  to  a  geocentric  system.  The  new  coordinate  system  uses  the  axes  £, 
r|,  and  where  the  ^-axis  is  pointing  towards  the  first  observed  position,  the  direction  of 


26 


the  third  observation  intersects  the  positive  rj-axis,  and  C,  is  mutually  perpendicular  to 
both  %  and  r|  as  seen  in  Figure  9. 


11 


Fig.  9.  —  This  illustrates  the  coordinate  transformation. 


In  this  case, 


Pix 


(p3XPl) 


|pi  xG 

p3(pi 

Jjxp,) 

•Pl)"Pl(l 

V 

a  a 

P3"P 

-  (P3  *  P. 

l(p3#Pl) 

V1  “(P3  *Pl)2 


and 

/\ 

A  c#  A 

9=^xn 


(18) 


The  next  step  is  making  a  rotation  matrix  that  will  change  the  axis  systems.  The  unit 


vectors  in  the  direction  of  the  observation  become, 
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(19) 
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A  first  approximation  to  the  constants  cj  and  C3  is  found  by  only  taking  the  first  terms 
in  the  f  and  g  series 


f  =  1  and  g  =  t;.  - 12 


and  by  allowing 

T,  =(t,  -t2),  T  2  =  (t3  -  tj),  T3  =(t3-t2) 
Substituting  these  values  into  equation  (15)  gives 


(20) 


(21) 


c  _ _ l3  l2 _ _  l3 

l(t3  -t2)-(ti  -^2)! 

and  (22) 

c  hzh  -  Ti 

l(t3-t2)-(t,-t2)l  T3-V 

Next,  the  pi ,  p2,  and  p3  components  of  equation  (16)  are  derived  from  equation  (19), 
which  yields 


-C]Zi  +Z2  -C3Z3 


_  p2p2  +c,Y|  -  Y2  +C3Y3 

P3  “  ’ 

c3|X3 

_  P2^-2  ~C3P3^3  ~  -^2  ~*"C3-^3 


(23) 


c, 

where  X„  Y„  and  Z,  are  the  heliocentric  positions  of  the  Earth  at  time  i.  These  three 
equations  with  three  unknowns  are  solved  producing  pi,  p2,  and  P3.  Once  the  magnitudes 
are  found,  r,  can  be  found  by  using  equation  (15).  But  ci  and  C3  also  depend  on  r,-, 


|r2  x  r3 

|ri  x  r3 

ri  X  I3 


(24) 
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These  are  Gauss'  sector  triangle  ratios.  The  constants,  Ci  and  C3,  are  now  solved 
iteratively  until  the  values  for  r(t)  converge.  Finally,  r  can  be  found  by  using  equation 
(17). 


5.2  Determination  of  Orbital  Element^ 

Next,  the  values  of  r  and  rare  used  to  find  the  orbital  elements  (e.g.,  Danby 
1995).  Since  Kepler's  First  Law  states  that  all  orbits  are  ellipses,  the  orbital  elements  are 
values  used  to  describe  an  ellipse.  Therefore,  a  combination  of  the  properties  of  an 
ellipse,  conservation  of  energy,  and  Kepler's  Laws  are  used  to  determine  the  orbital 
elements.  To  find  the  semi-major  axis  (a)  the  equation 


a  r  p 

is  used.  Next,  to  find  the  eccentricity  (e)  two  equations  with  two  unknowns  can  easily  be 
solved  for, 


r  «r 


e  cos(e)  =  1 —  and  e  sin(E)  =  — 
a  y  ap 


The  eccentric  anomaly  (E)  is  defined  in  Figure  10. 


(26) 


Fig.  10.  --  Illustration  of  the  definition  of  the  eccentric  anomaly  (E).  Every  ellipse  has  an  associated  circle. 
In  this  figure,  p  is  the  corresponding  point  on  the  circle  to  the  arbitrary  point  P  on  the  ellipse. 
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Now  that  both  e  and  E  are  known,  the  mean  anomaly  (Mo)  can  be  found  using 

M0=E-esin(E)  (27) 


The  longitude  of  the  ascending  node,  Q,  the  inclination  angle,  i,  and  the  argument  of 
perihelion,  co,  can  be  found  from  the  x,  y,  and  z  components  of  h, 

hx  =hsin(Q)sin(z), 

hy  =hcos(Q)sin(z),  (28) 

hz  =  h  cos(z) 

where  h  is  the  angular  momentum  per  unit  mass,  h  =  r  x  r ,  and  is  perpendicular  to  the 
orbital  plane.  Finally,  the  argument  of  perihelion  (co)  needs  to  be  found.  The  vector 
pointing  towards  the  perihelion  with  an  absolute  value  e  is 

e  =  (ex,e  ,e  )  = - r  (29) 

E 

Also,  from  geometry, 

e  cos(co)  =  ex  cos(q)+  ey  sin(Q)  (30) 


and 


esin(co)=-7-y^ 

sin(z') 


(31) 


By  combining  these  three  equations,  the  argument  of  perihelion  (co)  can  be  found. 


3.3  Perturbed  Motion 
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Gauss’  Method  only  takes  the  gravitational  interaction  with  the  Sun  into  account. 
However,  other  objects  in  the  Solar  System,  such  as  Earth,  may  influence  the  orbit  of  the 
NEA.  Perturbing  forces  are  the  gravitational  forces  exerted  on  the  asteroid  by  bodies 
other  than  the  Sun.  These  forces  cause  the  orbit  to  change  over  time.  There  are  two  main 
methods  that  are  used  to  account  for  these  forces  explicitly.  The  first  is  based  on  tidal 
perturbations  in  a  binary  system  (such  as  the  Earth-Moon  system).  It  was  first  introduced 
by  Newton  and  formalized  by  Lagrange  as  partial  derivatives  of  the  perturbing  force  with 
respect  to  the  elements.  The  second  method  comes  from  Gauss.  Gauss  resolves  the 
perturbing  force  into  a  coordinate  system,  which  is  rotating  with  the  orbital  motion  (e.g., 
Danby  1988).  However,  these  calculations  are  complicated,  computer  intensive  and 
beyond  the  scope  of  this  project,  especially  considering  that  it  is  fruitless  to  add 
perturbations  to  an  orbit  if  more  observations  are  needed  to  solidify  the  values  of  the 
orbital  elements.  On  the  other  hand,  the  Minor  Planet  Center  calculates  orbits  including 
perturbations  where  it  is  justified.  The  conclusions  made  here  regarding  the  possibility  of 
a  close  encounter22  with  these  NEAs  do  not  differ  significantly  from  those  made  by  the 
Minor  Planet  Center. 


22 


A  close  encounter  is  within  0.2  AU  of  Earth. 


4.  Lightcurves 
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A  lightcurve  is  obtained  by  observing  a  source  regularly  throughout  at  least  one 
foil  night  from  sunset  to  sunrise.  As  shown  in  Figure  25,  a  lightcurve  is  a  graph  of 
magnitude  vs.  time  used  to  determine  three  physical  parameters.  First,  a  lightcurve  can 
be  used  to  determine  the  general  shape  of  the  asteroid.  As  the  asteroid  rotates,  different 
sides  face  the  Earth.  Assuming  the  asteroid  is  not  a  sphere  and  that  its  rotational  axis  is 
not  along  our  line  of  sight,  as  it  rotates  the  reflected  light  from  the  Sun  will  vary.  So  by 
looking  at  the  amplitude  of  magnitude  variations,  the  shape  of  the  asteroid  is  estimated. 

W 

The  flux  of  light  observed  at  the  Earth  (in  — -)  is  given  by 

m 


f  =  L„  — — As  7t 


1  .  2  1 


4r  7t  4p  7i 


(32) 


where  L.  is  the  luminosity  of  the  Sun  (in  W),  A  is  the  albedo23  of  the  asteroid,  r  is  the 
distance  from  the  Sun  to  the  asteroid,  p  is  the  distance  from  the  Earth  to  the  asteroid,  and 
s  is  the  radius  of  the  asteroid.  As  the  asteroid  rotates,  the  cross-sectional  area  facing  the 
Earth  changes,  but  all  other  parameters  remain  the  same24.  So  the  maximum  and 
minimum  flux  measured  is  related  to  the  change  in  magnitude  as 


23  Albedo  is  the  fraction  of  incident  light  that  is  reflected. 

24  All  terms  in  the  two  fluxes  are  common  terms  except  for  the  asteroid’s  radius.  The  luminosity  of  the  Sun 
is  a  constant.  The  distance  of  Earth  from  the  Sun  and  the  distance  of  the  asteroid  from  the  Sun  can  both  be 
assumed  to  be  constant  since  only  a  small  amount  of  time  passes  from  one  observation  to  the  next. 
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mmax  ‘  mmin  =  T^g 


f 

Vxmax  J 


=fi°gi 


( s  .  v 

amin 

V^niax  J 


(33) 


(0  \ 


=  5  log 


V^max  J 


Then  solving  for  the  ratio  of  the  radii  gives  the  relative  shape  of  the  asteroid, 


[m  max  ~m  min  1 

5  ■  (34) 

o 

max 

A  lightcurve  is  also  used  to  determine  the  rotational  period  of  the  asteroid, 
assuming  that  every  other  peak  in  the  lightcurve  is  due  to  reflected  sunlight  off  the  same 
face  of  the  asteroid.  The  rotational  period  is  obtained  by  finding  a  pattern  in  the 
lightcurve  that  repeats  several  times.  The  time  between  the  two  patterns  is  the  rotational 
period. 

Lightcurves  of  an  asteroid  from  at  least  three  different  times  that  the  asteroid  is  at 
opposition25  are  needed  in  order  to  find  an  asteroid’s  pole  axis  (e.g.,  Kiss,  Szabo,  & 
Sameczky  1999).  Since  this  would  take  several  years,  this  project  focused  on  finding  the 
rotational  period  and  general  shape  of  an  asteroid  from  its  lightcurve  without  finding  its 
pole  axis.  Both  the  rotational  period  and  shape  are  important  parameters  if  a  plan  is 
required  to  deflect  the  asteroid  from  impacting  Earth. 


25  Opposition  is  the  point  in  the  object's  orbit  when  Earth  is  between  it  and  the  Sun,  making  it  the  brightest 
that  we  see  it. 


5.  Results  and  Discussion 
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5. 1  Cooperation  Process 

Ideally,  an  asteroid  should  be  observed  on  nearly  a  dozen  nights  spread  out  over 
three  months  in  order  to  calculate  an  orbit  reliable  for  decades.  For  a  number  of  reasons, 
it  is  not  possible  for  a  single  observatory  to  obtain  enough  data  to  do  this  alone. 

A  major  problem  almost  every  optical  observatory  faces  is  bad  weather,  such  as 
rain,  snow,  or  cloud  cover.  This  severe  weather  obviously  closes  the  observatory  for  the 
night,  but  even  more  mild  conditions  threaten  observing.  For  example,  a  light  haze  will 
impede  observing  because  the  particles  in  the  clouds  scatter  the  light  from  astronomical 
sources.  The  asteroids  are  so  faint  that  any  blurring  caused  by  scattered  light  makes  them 
disappear  into  the  ambient  noise  of  the  image.  Another  hazardous  weather  condition  is 
the  humidity.  If  the  humidity  is  too  high,  then  condensation  can  form  on  the  CCD. 
Condensation  on  the  CCD  can  cause  a  couple  of  problems.  First,  it  scatters  the  light 
hitting  the  CCD,  causing  the  image  to  be  blurred  and  faint  if  at  all  visible.  Also,  if 
condensation  starts  to  form  in  the  control  box,  then  condensation  on  the  wires  can  cause 
short  circuits.  Of  course,  even  with  ideal  weather  conditions,  equipment  failure  is  still  a 
threat.  During  the  course  of  this  project,  the  CCD  broke  twice,  the  computers  in  the 
dome  were  replaced,  and  the  telescope  tracking  failed.  Also,  if  the  asteroid  happens  to  be 
moving  away  from  Earth  while  it  is  being  tracked,  then  it  will  get  fainter.  This  may  make 
it  impossible  to  track  the  asteroid  for  the  full  three  month  period. 

One  of  the  reasons  that  the  Minor  Planet  Center  exists  is  because  no  single 
observatory  can  reliably  obtain  enough  observations  of  an  asteroid  to  determine  the  orbit. 


Therefore,  whenever  an  observatory  finds  the  position  of  an  NEA,  it  reports  this 
position  into  the  Minor  Planet  Center.  The  Minor  Planet  Center  then  posts  these 
coordinates  on  the  world-wide  web  so  that  anyone  may  use  them  to  calculate  the  best 
orbit  possible.  To  illustrate  why  it  is  important  to  use  as  much  data  as  possible.  Table 
gives  three  different  sets  of  orbital  elements  for  asteroid  2000  EB107.  One  orbit  is 
calculated  from  data  obtained  solely  at  the  U.S.  Naval  Academy,  the  second  orbit  is 
calculated  only  from  data  posted  on  the  world-wide  web  by  the  Minor  Planet  Center, 
while  the  third  orbit  is  calculated  from  both  of  these  data  sets. 


TABLE  1 

Orbital  Elements  for  2000  EB107 


Source 

a 

(AU) 

e 

(deg) 

Q 

(deg) 

(0 

(deg) 

M0 

USNA 

1.2 

0.63 

37 

180 

230 

64 

MPC 

3.1 

0.59 

26 

180 

43 

350 

USNA+MPC 

2.7 

0.55 

24 

180 

44 

340 

Fig.  11.  —  Bar  Graph  comparing  the  orbital  elements  of  2000  EB107. 
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As  can  be  seen  in  Figure  11,  there  are  great  differences  in  the  orbital  elements, 
especially  for  the  semi-major  axis  and  the  argument  of  perihelion.  The  semi-major  axes 
differ  by  more  than  1  AU,  while  the  difference  in  the  argument  of  perihelion  is  almost 
200°.  In  Figure  12,  these  three  orbits  are  shown.  Even  a  quick  inspection  of  Figure  12 
clearly  reveals  that  these  differences  in  the  orbital  elements  drastically  change  the  shape 
of  the  orbit. 


North 


Fig.  12.  --  The  orbit  of  2000  EB107  obtained  by  using  data  from  different  sources. 


Because  it  uses  the  largest  range  of  dates  over  which  the  asteroid  was  observed,  the  best 
orbit  came  from  the  combination  of  USNA  and  MPC  data  sets. 


5.2  Using  an  Astronomical  Calibrator^ 
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The  photometry  software,  Fitsblink,  does  not  necessarily  find  the  correct 
independent  magnitude  of  the  objects  in  an  image.  However,  it  does  give  the  correct 
relative  magnitude,  meaning  that  the  difference  in  magnitude  of  two  objects  in  an  image 
is  correct.  However,  the  difference  in  magnitude  of  an  asteroid  that  is  found  between  two 
images  is  not  only  due  to  an  actual  change  in  its  magnitude,  but  also  is  due  to  an  absolute 
calibration  error  between  images.  In  order  to  correct  for  this,  a  reference  star  is  used. 

This  is  simply  a  star  that  is  found  in  the  majority,  if  not  all,  of  the  images.  In  principle, 
the  star’s  magnitude  in  each  of  these  images  should  be  the  same,  but  it  is  not.  To 
calibrate  a  series  of  images,  all  magnitudes  are  referred  to  an  arbitrary  observation  of  the 
reference  star  in  one  of  the  images.  Figure  13  shows  the  lightcurve  of  2000  EE104  with 
uncalibrated  magnitudes  and  calibrated.  In  particular,  points  near  6  hours  UT  have 
changed  by  three  magnitudes.  This  calibration  greatly  influences  the  interpretation  of  the 
data.  As  a  final  note,  the  main  source  of  error  in  the  magnitudes  is  due  to  the  Guide  Star 
Catalog.  The  Guide  Star  Catalog  only  reports  magnitudes  to  within  one  third  of  a 
magnitude.  This  large  error  effectively  drowns  out  any  other  source  of  error  so  that  only 
this  error  need  be  addressed. 


26  UT  is  the  universal  time  or  the  local  standard  time  at  the  Prime  Meridian. 


Fig.  13.  —  Comparison  of  a  raw  and  calibrated  lightcurve. 
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5.5  Orbital  Results  and  Discussion 

For  this  project,  once  enough  observations  were  made  to  obtain  an  orbit,  the 
orbital  elements  were  found  using  a  program  provided  by  Spacewatch.  This  program 
uses  five  different  methods  to  find  the  orbital  elements  and  computes  the  RMS  of  the 
fits.  The  five  methods  are  Gauss’  (see  Section  3.1),  Circular,  and  three  Viasala 
calculations.  The  Circular  and  three  Viasala  methods  calculate  the  orbital  elements 
assuming  the  eccentricity  to  be  0,  0.15,  0.30,  and  0.45  respectively.  Circular  or  one  of  the 
three  Viasala  calculations  is  used  when  there  is  very  little  data,  so  the  eccentricity  must 
be  assumed,  and  only  the  asteroid’s  position  for  the  following  night  is  required.  The  best 
fitting  orbit  has  the  lowest  RMS  value.  Ideally,  Gauss’  Method  is  used  because  it  does 
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not  assume  any  orbital  elements.  Table  2  displays  the  orbital  elements  that  were  found 


for  each  asteroid. 


TABLE  2 

Orbital  Elements  for  the  Asteroids  that  were  Tracked 


Asteroid 

Calc. 

Method 

a 

(AU) 

e 

(deg) 

Q 

(deg) 

0) 

(deg) 

M0 

Range  of 
Dates 

2000 

DM8 

Gauss 

1.2 

0.41 

18 

340 

300 

260 

7 

2000 

DH8 

Gauss 

2.4 

0.45 

42 

160 

21 

350 

33 

2000 

DOl 

Viasala 

2.4 

0.45 

120 

340 

180 

360 

1 

2000 

EB107 

Gauss 

2.7 

0.55 

24 

180 

44 

350 

24 

2000 

DQ110 

Gauss 

3.3 

0.63 

58 

350 

95 

20 

13 

2000 

ED  104 

Gauss 

1.3 

0.28 

31 

190 

230 

110 

23 

2000 

EE14 

Viasala 

1.3 

0.15 

98 

180 

43 

wm 

27 

2000 

EB14 

Viasala 

1.5 

0.15 

99 

160 

■fl 

4 

2000 

EDM 

Gauss 

0.87 

0.23 

4.4 

5.5 

310 

4 

2000 

DP  107 

Gauss 

1.4 

0.37 

8.4 

360 

290 

260 

31 

2000 

EE  104 

Gauss 

0.98 

0.18 

3.1 

27 

280 

260 

22 

Note-As  discussed  in  the  text  there  are  two  sources  of  systematic  error:  gravitational  perturbations  and 
insufficient  observations  of  the  orbit.  The  random  errors  in  the  astrometry  are  approximately  ±1  arcsec  and 


±5  s. 


27  RMS  is  the  root  mean  squared.  A  low  root  mean  squared  means  the  data  lie  close  to  the  orbital  fit  and  is 
therefore  desirable. 
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Asteroid  2000  DM8 

Gauss’  Method  was  used  to  determine  the  orbital  elements  for  2000  DM8.  From 
the  orbit  it  appears  to  be  an  Apollo  asteroid.  As  can  be  seen  in  Figure  14,  this  initial 
calculation  does  not  predict  a  close  approach  to  Earth.  At  the  places  where  it  appears  to 
cross  Earth’s  orbit,  it  is  far  above  Earth’s  orbit  due  to  its  inclination  angle  of  18.3°. 

However,  these  conclusions  are  preliminary.  Asteroid  2000  DM8  was  lost  after 
seven  days  of  tracking  because  it  became  too  faint  for  the  equipment  to  detect. 


Fig.  14.  -  Orbit  of  2000  DM8. 
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Asteroid  2000  DH8 

The  orbital  elements  of  2000  DH8  were  also  found  using  Gauss'  Method.  As  can 
be  seen  in  Figure  15,  it  does  not  appear  to  come  close  to  the  orbit  of  Earth.  This  orbit  is 
reliable,  since  it  was  based  on  more  than  a  month’s  worth  of  observations.  Because  this 
asteroid  resides  primarily  in  the  main  asteroid  belt28  and  is  still  bright  enough  to  be  easily 
seen,  its  size  must  be  at  least  on  the  order  of  50  km.  Therefore,  it  should  continue  to  be 
tracked  through  May  so  that  its  orbit  may  be  determined  more  accurately. 

Fig.  15. -Orbit  of  2000  DH8. 


28  The  main  asteroid  belt  is  between  the  orbits  of  Jupiter  and  Mars  and  is  where  most  asteroids  can  be 
found. 
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Asteroid  2000  DPI 

Figure  16  displays  2000  DOl’s  orbit.  This  orbit  was  based  only  on  one  night's 
observation.  Therefore,  a  Viasala  calculation  that  assumed  an  eccentricity  of  0.45 
produced  the  best  orbit.  This  orbit  can  only  give  the  approximate  position  of  the  asteroid 
for  a  few  nights;  it  cannot  be  used  to  determine  if  the  asteroid  may  have  a  close  approach 
with  Earth.  For  these  reasons,  more  observations  need  to  be  made  in  order  to  obtain  a 
Gaussian  orbit.  However,  it  is  likely  to  be  lost  because  soon  it  will  be  too  close  to  the 
Sun  to  observe. 

Fig.  16. -Orbit of 2000 DOl. 
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Asteroid  2000  EB107 

A  fairly  well  determined  orbit  for  2000  EB107  is  displayed  in  Figure  17.  It 
appears  to  be  an  Amor  asteroid.  From  the  figure,  it  seems  that  the  asteroid  will  not 
collide  with  the  Earth  due  to  the  fact  that  at  the  only  place  that  its  orbit  comes  within 
Mars'  orbit,  it  is  much  farther  out  of  the  ecliptic  with  an  inclination  angle  of  24°. 
Therefore,  this  asteroid  is  not  a  high  priority  for  follow-up  observations. 

Fig.  17.  -  Orbit  of  2000  EB107. 


North 
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Asteroid  2000  DQ1 10 

Despite  the  fact  that  the  orbit  of  2000  DQ1 10  was  found  by  Gauss’  Method,  the 
most  reliable  method  used,  it  is  still  suspect.  Like  2000  DH8,  this  asteroid  appears  to  be 
a  main  belt  asteroid,  rather  than  an  NEA  (Figure  1 8).  It  has  a  large  semi-major  axis  and 
eccentricity.  Also,  it  was  only  observed  over  a  period  of  thirteen  days,  requiring  that 
further  tracking  be  done  to  obtain  a  trustworthy  orbit. 

Fig.  18.  -  Orbit  of  2000  DQ1 10. 


North 


Asteroid  2000  ED  104 


With  data  spanning  twenty-four  days,  a  decent  orbit  of  2000  EDI 04  was 
determined  using  Gauss’  Method.  Figure  19  suggests  that  2000  ED  104  is  an  Apollo 
asteroid  and  that  it  might  have  a  close-approach  with  Earth.  In  fact,  the  Minor  Planet 
Center  predicts  that  it  will  be  within  0.17AU  of  Earth  on  9  October  2002. 


Fig.  19.  --  Orbit  of  2000  ED104. 


Asteroid  2000  EE14  and  2000  EB14 
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Although  2000  EE14  was  observed  over  a  time  of  twenty-seven  days,  a  Gaussian 
solution  could  not  be  calculated  due  to  the  size  of  the  gaps  between  observations. 
Asteroid  2000  EB14  was  tracked  for  only  four  days.  For  both  of  these  asteroids,  the  fact 
that  a  Viasala  solution  had  to  be  used  and  their  high  inclination  angles,  as  seen  in  Figures 
20  and  21,  make  these  asteroid's  orbit  untrustworthy.  Both  asteroids  require  more 
observations  so  that  a  Gaussian  orbit  can  be  calculated. 

Fig.  20.  -  Orbit  of  2000  EE14. 
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Asteroid  2000  ED  14 

Although  2000  ED  14  has  only  been  observed  for  a  period  of  four  days,  its 
Gaussian  orbit  is  tantalizing.  Both  the  orbit  shown  in  Figure  22  and  the  Minor  Planet 
Center  agree  that  it  may  have  a  close  approach  with  Earth.  In  fact,  according  to  the 
Minor  Planet  Center,  this  asteroid  was  merely  0.04AU  from  Earth,  twice  the  Moon’s 
distance  from  Earth,  on  March  16, 2000  (spring  break),  and  it  is  predicted  to  come  within 
0.2AU  of  Earth  in  March  2013,  March  2016,  and  March  2029.  Clearly  this  asteroid  is  a 
top  priority  for  further  investigation. 

Fig.  22.  -  Orbit  of  2000  EDM. 

North 


Neptune 
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Asteroid  2000  DP  107 

This  fairly  well  determined  orbit,  displayed  in  Figure  23,  shows  a  close  encounter 
with  Earth.  However,  the  Minor  Planet  Center  does  not  report  that  this  asteroid  has  a 
close  approach  with  Earth.  Due  to  this  discrepancy  and  the  fact  that  it  may  be  a  threat  to 
Earth,  further  observations  of  2000  DP  107  is  a  top  priority. 


Fig.  23.  -  Orbit  of  2000  DP  107. 


Asteroid  2000  EE104 
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Both  Figure  24  and  the  Minor  Planet  Center  agree  that  this  orbit  will  have  a  close 
approach  with  Earth.  In  fact,  it  was  calculated  that  only  six  days  before  the  due  date  of 
this  project  (1 1  April  2000)  this  asteroid  was  0.04AU  from  Earth.  It  is  predicted  that  this 
asteroid  will  have  a  close  encounter  with  Earth  twelve  times  in  the  next  thirty-three  years. 
Clearly  this  asteroid  is  a  top  priority  for  further  investigation. 


Fig.  24.  ~  Orbit  of  2000  EE  104. 


5.4  Lightcurve  Result s  and  Discussion 
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This  section  presents  the  lightcurves  that  were  found  for  five  of  the  eleven 
asteroids  that  were  being  tracked.  Each  lightcurve  was  observed  over  the  period  of  one 
full  night  from  sunset  to  sunrise.  Two  asteroids  were  observed  for  an  additional  night. 
On  1  April  2000  lightcurves  were  obtained  for  asteroids  2000  EB107  and  2000  EE  104, 
while  2000  DM8  and  2000  DH8  were  observed  on  4  March  and  6  March  2000.  A  partial 
lightcurve  was  observed  for  2000  EE14  on  30  March  2000.  Each  image  was  then 
calibrated  and  the  lightcurves  obtained. 


Asteroids  2000  EB107  and  2000  EE  104 

The  lightcurves  for  2000  EB107  and  2000  EE  104  are  displayed  in  Figure  25. 
Except  for  a  gap  of  two  hours  during  which  time  other  asteroids  were  being  tracked, 
observations  were  made  on  each  asteroid  every  twenty  minutes  from  when  the  asteroid 
first  rose  to  when  it  set.  Due  to  these  gaps  in  the  data,  only  limits  for  the  period  can  be 
found.  The  time  from  one  maximum  to  the  next  is  at  most  half  of  the  period.  Therefore 
the  lower  limits  of  the  rotational  period  are  3.6  hours  for  2000  EB107  and  8  hours  for 
2000  EE  104 

For  2000  EB107,  the  brightest  magnitude  is  17.6  and  faintest  is  19.5.  The 
magnitude  range  for  2000  EE104  is  from  15.2  to  16.1.  Using  equation  (34), =  2.4 

Smin 

for  2000  EB107  and  =  1 .5  for  2000  EE104.  These  ratios  are  illustrated  in  Figure  27. 

®min 


Asteroids  2000  DM8  and  2000  DH8 
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The  data  for  these  lightcurves,  Figure  26,  were  taken  over  two  nights.  The 
rotational  period  of  2000  DM8  is  4.6  hours.  Because  a  clear  pattern  is  not  established  in 
2000  DH8  the  rotational  period  can  only  be  pinned  down  to  a  range  between  3.2  and  5.5 
hours.  The  brightest  and  faintest  magnitudes  for  2000  DM8  are  17.1  and  19.2.  This 

gives  =  2.6 .  The  brightest  and  faintest  magnitudes  for  2000  DH8  are  17.6  and  19.0. 

Smin 

g 

This  gives  -2152-  =  1.9.  These  ratios  are  represented  by  the  shapes  seen  in  Figure  27. 

^min 

Only  about  one  and  a  half  nights  worth  of  data  was  used  for  these  lightcurves  due  to  the 
fact  that  the  tracking  on  the  telescope  began  to  malfunction  in  the  middle  of  the  second 
night.  Otherwise,  these  asteroids  were  observed  every  twenty  minutes  on  both  nights. 
Since  the  time  of  night  on  both  nights  do  not  correspond  to  the  same  phase  of  rotation  of 
the  asteroid,  the  data  from  the  second  night  was  shifted  on  the  x-axis  so  that  it  best  fit  the 
first  night’s  data. 
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Fig.  27.  --  Rough  shapes  based  off  radius  ratios  for  2000  EB107,  2000  EE104, 2000  DH8,  and 

2000  DM8. 


2000  EB107  2000  DH8 


2000  EE 104  2000  DM8 


Asteroid  2000  EE  14 

The  lightcurve  for  2000  EE14  is  shown  in  Figure  28.  As  can  be  seen  from  the 
plot,  nothing  conclusive  can  be  determined.  High  humidity  during  the  night  caused 
condensation  to  form  on  the  window  of  the  CCD.  Before  this  problem  was  discovered, 
2000  EE  14  was  observed  for  approximately  three  hours.  Most  likely  these  data  are 
corrupted  by  condensation.  This  plot  is  included  to  illustrate  some  of  the  observational 
problems  that  were  encountered  during  this  project  and  in  hopes  of  helping  others  who 
may  do  similar  observations  at  USNA. 


6.  Conclusion 
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This  project  has  emphasized  the  importance  of  carefully  tracking  NEAs.  Of  the 
eleven  asteroids  that  were  tracked  for  this  project,  approximately  half  will  make  a  close 
approach  to  the  Earth.  Clearly  these  asteroids  require  further  monitoring.  But  more 
importantly,  the  estimated  1500  undiscovered  NEAs  will  need  careful  tracking. 

Small  observatories  like  the  now  fully  operational  one  on  the  roof  of  Michelson 
Hall  are  ideally  suited  for  such  work.  Many  lessons  have  been  learned  as  a  result  of  this 
project.  First,  observing  in  the  summer  in  Annapolis  is  nearly  impossible  due  to  the  high 
humidity.  Second,  it  is  not  useful  to  take  lightcurves  and  track  asteroids  on  the  same 
night.  This  dual  procedure  leads  to  gaps  in  the  lightcurve  data  which  limit  the  analysis. 
Third,  in  order  to  compensate  for  degradation  in  tracking,  the  telescope  must  be  re¬ 
initialized  every  one  to  two  hours  so  that  the  desired  object  does  not  leave  the  CCD's 
field  of  view.  Finally,  same  night  analysis  of  observations  must  be  made.  This  is  so  that 
any  errors  in  the  data  can  be  seen  early  on  and  corrected,  such  as  tracking,  condensation. 


or  focus. 


7.  Future  Work 
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For  obvious  reasons,  this  research  is  an  ongoing  project  requiring  more  work. 
First,  the  top  priority  asteroids,  those  with  close  approaches  to  Earth,  will  continue  to  be 
tracked  so  that  a  more  definite  orbit  may  be  calculated.  In  addition,  since  the  number  of 
unknown  NEA  orbits  is  still  large,  observations  of  new  asteroids  will  be  done  so  that  an 
orbit  may  be  calculated  for  them.  Not  only  will  I  do  this  over  the  next  few  months,  but 
MIDN  Aaron  Fielder,  USN  will  be  doing  this  as  part  of  an  independent  research  project 
next  semester.  Also,  more  lightcurve  data  will  be  taken  of  2000  EB107  and  2000  EE  104 
in  order  to  fill  in  the  existing  gaps  in  their  lightcurves.  This  is  important  because  if 
MIDN  Fielder  then  observed  these  same  asteroids  next  year  at  opposition,  only  one  more 
lightcurve  would  be  needed  to  determine  the  pole  axes  of  the  asteroids.  Another  thing 
that  will  be  done  over  the  summer  is  a  consolidation  of  the  telescope  workstation  to  a 
single,  remote  observing  site  on  the  third  floor  of  Michelson  Hall.  Currently,  the 
computer  which  controls  the  telescope  is  in  one  room,  the  computer  which  controls  the 
CCD  is  on  the  roof  of  Michelson  Hall,  and  the  computer  that  analyzes  the  data  is  in  a 
third  room.  The  main  obstacle  to  this  consolidation  will  be  finding  a  way  to  control  the 
dome  remotely  and  accurately.  A  remote  site  is  important  because  the  dust  and  weather 
conditions  in  the  dome  have  already  broken  one  monitor  and  keyboard. 
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Appendix  A:  Laplace’s  Method 
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Laplace's  Method  for  determining  orbital  elements  uses  four  equations.  Three  of 
these  equations  are  derived  from  the  asteroid’s  geometrical  relationship  with  Earth  and 
the  Sun,  while  the  fourth  represents  the  gravitational  attraction  between  the  Sun  and  the 
asteroid.  The  geometrical  relationship  between  the  asteroid,  Earth,  and  the  Sun  can  be 
seen  in  Figure  8.  Summing  the  vectors,  gives  equation  (15) 

R  +  r  =  p  =  pp.  (15) 

Taking  the  time  derivative  of  equation  (15)  yields, 

R  +  r  =  pp  +  pp.  (35) 

Taking  the  time  derivative  once  more  gives, 

R  +  r  =  pp  +  pp  +  pp  +  pp  ^ 

=  pp  +  2pp  +  pp 

The  third  equation  relates  the  gravitational  force  between  the  bodies.  Both  the  Earth  and 
the  asteroid  have  a  gravitational  attraction  to  the  Sun,  since  the  mass  of  the  Sun  is  so 
much  larger  than  the  masses  of  the  Earth  or  of  the  asteroid,  the  masses  of  the  Earth  and 
the  asteroid  are  ignored.  This  relationship  gives 

r  =  -p-^-  and  R  =  -p— r  (37) 

r3  R3 

where  p=GMsun.  If  these  two  equations  are  summed  and  equation  (15)  is  substituted  for  r 


then. 


By  using  the  Law  of  Cosines, 

r2  =  p2  +  R2  -  2Rp  cos(\]/)  (40) 

is  obtained.  Therefore,  the  four  equations  needed  for  determining  the  orbital  elements 
using  Laplace’s  Method  are  equations  (15),  (35),  (39),  and  (40). 

Three  observations  of  the  asteroid  give  three  values  for  p  and  t.  From  these 

observations,  approximate  values  for  p  and  p  are  estimated.  The  initial  approximation 
used  in  Laplace's  method  involves  finding  values  for  p  and  p  at  some  point  in  time.  If 
this  point  in  time  is  assumed  to  be  to,  then  for  small  t,  p  can  be  expanded  as  a  Taylor 
series, 

P  =  (p),.„+t(p),+L!  (£)+...  (41) 

By  stopping  this  series  after  the  third  term,  only  three  observations  are  necessary  to  find 

( p  )0  and  ( p  )o.  It  is  best  to  take  t0  as  either  the  arithmetic  mean  of  the  times  or,  if  there 
are  only  three  observations,  as  the  central  time. 


The  next  two  steps  are  taken  in  order  to  reduce  the  number  of  unknown 
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variables.  First,  the  dot  product  of  equation  (39)  is  taken  with  (p  x  p),  so  that 


Next,  the  dot  product  of  equation  (39)  is  taken  with  (p  x  p), 


(42) 


R3 


(pxp).R  =  -2p(pxp).p, 

) 


(43) 


Once  this  is  done,  equation  (40)  and  equation  (42)  can  be  solved  for  r  and  p,  equation 


(43)  can  then  be  solved  for  p ,  and  equation  (15)  and  equation  (35)  can  be  solved  for  r 

90 

and  r ,  from  which  the  elements  of  the  orbit  can  be  calculated. 


29  Solving  for  the  final  values  of  r  and  r  is  a  complicated  process  involving  a  transcendental  equation. 
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